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Abstract. Using modern nucleon-nucleon interactions in the description of the A = 3, 4 nuclei, it is not possible 
to reproduce both the three- and four-nucleon binding energies simultaneously. This is one manifestation of the 
necessity of including a three-nucleon force in the nuclear Hamiltonian. Several models of the three-nucleon 
force exist and are applied in the description of light nuclei. However, as it is discussed here, a simultaneous 
description of the three- and four-body binding energies and the n - d doublet scattering length seems to be 
problematic. Accordingly, a comparative study of some of these models is performed. In a different analysis, we 
study applications of the Kohn Variational Principle, formulated in terms of integral relations, to describe N — d 
scattering processes. 



1 Introduction 

Realistic nucleon-nucleon (NN) potentials reproduce the 
experimental NN scattering data up to energies of 350 Me V 
with a x 1 per datum close to 1 . However, the use of these 
potentials in the description of the three- and four-nucleon 
bound and scattering states gives a x 1 P er datum much 
larger than 1 (see for example Refs. 01121 ). In order to im- 
prove that situation, different three-nucleon force (TNF) 
models have been introduced so far. Widely used in the 
literature are the Tucson-Melbourne (TM) and the Urbana 
IX (URIX) models [3 4|. These models are based on the 
exchange mechanism of two pions between three nucle- 
ons with the intermediate excitation of a A resonance. The 
TM model has been revisited within a chiral symmetry ap- 
proach [5 1, and it has been demonstrated that the contact 
term present in it should be dropped. This new TM poten- 
tial, known as TM', has been subsequently readjusted [6|. 
The final operator structure coincides with that one given 
in the TNF of Brazil already derived many years ago Q. 
Recently, TNFs have been derived based on chiral effective 
field theory at next-to-next-to-leading order (8]. A local 
version of these interactions (hereafter referred as N2LOL) 
can be found in Ref. [9 1. All these models contain a certain 
number of parameters that fix the strength of the interac- 
tion. It is a common practice to determine these parameters 
from the three- and four-nucleon binding energies. A par- 
ticular TNF is in general associated to a specific NN poten- 
tial and the sum of the two interactions forms the nuclear 
potential energy. The two- and three-nucleon interactions 
derived using chiral effective field theory are consistently 
constructed. However a particular TNF can be used asso- 
ciated with different NN interactions. As a consequence, 
the parametrization of a particular TNF could change since 
different NN potentials predict different A — 3,4 binding 
energies. 
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More recently, a new class of two-nucleon interactions 
has been obtained (Viow-k potentials). With the purpose of 
eliminating the high-momentum part of the interaction, the 
Hilbert space has been separated into low and high mo- 
mentum regions and the renormalization group method has 
been used to integrate out the high momentum components 
above a cutoff A IflOll . The value for A is typically chosen 
to reproduce the triton binding energy. 

All these potential models can be used to study bound 
and scattering states in the A = 3,4 systems in order to 
extract information about their capability to describe the 
nuclear dynamics. Besides the bound state energies, in the 
A — 3 system, the n - d doublet scattering length 2 a m { 
can give valuable information. In principle this quantity 
is correlated, to some extent, to the A = 3 binding en- 
ergy through the so-called Phillips line 111 11121 . However 
the presence of TNFs of the type studied here breaks this 
correlation. Therefore 2 a„d emerges as an independent ob- 
servable. Due to the lack of excited states in the A = 3 sys- 
tem, the zero energy state is the first one above the ground 
state. In the case of n - d scattering at zero energy, the 
J — \ state is orthogonal to the triton ground state and, 
for this reason, it presents a node in the relative distance be- 
tween the incident nucleon and the deuteron. The position 
of the node is related to the scattering length and it is also 
sensitive to the relation between the overall attraction and 
repulsion of the interaction. Several of the realistic NN po- 
tentials underestimate the triton binding energy. Therefore 
by adding a TNF, with the strength fixed for example to re- 
produce the triton binding energy, the balance between the 
total attraction and repulsion in the potential changes. This 
leads to a modification in 2 a„d and this modification de- 
pends on the parameters in the TNF. The determination of 
the TNF parametrization able to describe the triton binding 
energy B( 3 H), the a-particle binding energy B( 4 He) and 
2 a„d has been analyzed in Ref. [8] for a TNF derived from 
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chiral effective field theory. A similar analysis has not been 
done for the local TNF models URIX, TM' and N2LOL. 

In Refs. |[T3TT4l results for 5( 3 H), fi( 4 He), 2 a„ d are given 
using different combinations of NN interactions (see Ta- 
ble|T}. Those results indicate that the models are not able to 
describe simultaneously the A = 3,4 binding energies and 
2 a n d- In order to analyze further the mentioned discrepan- 
cies, here we study potential models constructed summing 
to the AVI 8 NN potential |[T5l the three-nucleon interac- 
tions of TM', URIX and N2LOL. Parametrizations of the 
URIX and TM' models already exist in conjunction with 
the AVI 8 potential. Conversely the N2LOL force has been 
constructed using the N3LO-Idaho potential from Entem et 
al. (T6'|. So, here we adapt its parametrization to reproduce, 
when summed to the AV 1 8 interaction, the triton binding 
energy. Different parametrizations of the three TNF mod- 
els are analyzed studying the description of fi( 3 H), fi( 4 He) 
and 2 a„d and some polarization observables in p-d scatter- 
ing. The calculations have been done using the hyperspher- 
ical harmonic (HH) method as given in Refs. 11171 1 81 1 9 20 1 
to describe bound and scattering states in A = 3,4 systems 
using local potentials. The extension to treat nonlocal po- 
tentials was given in Refs. H 1412 1 1 - 

In a different application devoted to study scattering 
states in few-nucleon systems, a discussion of the use of 
the integral relations derived in Ref. 11221 from the Kohn 
Variational principle (KVP) is given. It has been shown 
that starting from the KVP, the tangent of the phase-shift 
can be put in a form of a quotient where both, the nu- 
merator and the denominator, are given in the form of an 
integral relation. This is similar to what was proposed in 
Ref. [23 1, however its strict relation with the KVP has not 
been recognized. To be noticed that a general formulation 
of the scattering theory using surface-integrals is given in 
Ref. [24|. Here we would like to discuss some specific 
examples of the integral relations derived from the KVP. 
Starting the analysis in the simplest case, the A — 2 sys- 
tem, we show that they can be used to compute phase-shifts 
from bound state like functions. A second application of 
the integral relations regards the possibility of determin- 
ing phase-shifts from a calculation in which the Coulomb 
potential has been screened. All these examples serve to 
demonstrate the general validity of the KVP formulated in 
terms of integral relations. Due to their short-range nature, 
they are determined by the wave function in the interaction 
region and not from its explicit asymptotic behaviour. This 
means that each wave function solving {H-E)W — in the 
interaction region can be used to determine the correspond- 
ing scattering amplitude even if its asymptotic behaviour is 
not the physical one. 



2 The HH expansion for A = 3,4 systems 



In this section we briefly review the HH method for bound 
and scattering states. 



2.1 The HH Method for Bound States 

The nuclear wave function for the three-body system can 
be written as 

VP) = Yj C *\ W *) > (!) 

where | W^) is a suitable complete set of states, and p is an 
index denoting the set of quantum numbers necessary to 
completely specify the basis elements. 

The coefficients of the expansion can be calculated us- 
ing the Rayleigh-Ritz variational principle, which states 
that 

(d c ¥\H-E\¥) = 0, (2) 

where dcf indicates the variation of W for arbitrary in- 
finitesimal changes of the linear coefficients c M . Where the 
Hamiltonian of the system consists in the kinetic part plus 
two- and three-nucleon interaction terms 

i i< j i< j<k 

The problem of determining c M and the energy E is re- 
duced to a generalized eigenvalue problem, 

^<^|//-£|^><y =0. (4) 

The main difficulty of the method is to compute the ma- 
trix elements of the Hamiltonian H with respect to the ba- 
sis states yFfj). Usually H is given as a sum of terms (ki- 
netic energy, two-body potential, etc.). The calculation of 
the matrix elements of some parts of H can be more con- 
veniently performed in coordinate space, while for other 
parts it could be easier to work in momentum space. There- 
fore, it is important that the basis states | f^) have simple 
expressions in both spaces. The HH functions indeed have 
such a property. 

In the case of three nucleons of mass m the Jacobi vec- 
tors xi p ,X2 P correspond to a given particle permutation 
denoted with p, which specifies the particle order i, j, k, 

X2P = ^2 Vi ~ V>) ' 
[2 1 

x\ p = J -(r k - -(r; + rj)) . (5) 

Here p — 1 corresponds to the order 1,2,3. It is convenient 
to replace the modulii of X2 P and X\ p with the so-called 
hyperradius and hyperangle, defined as 

tan0 p = — '- . (7) 

X2p 

Note that p does not depend on the particle permutation 
p. The complete set of hyperspherical coordinates is then 
given by {p, £2^}, with 

= [xi p , x 2p ; <p p ] , (8) 
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and the suffix (p) recalls the use of the coordinate space. 
The expansion states | f^) of Eq. (Q]i are then given by 

\ 1 Pjf ) )=fi(pW{G}(^ ) ), (9) 

where fi(p) for I = 1 , . . . M is a complete set of hyperradial 
functions, chosen of the form 



fi(p) = r 



/! 



'a +5)! 



Lf\yp) e" 



(10) 



HereLfCyp) are Laguerre polynomials, and the non-linear 
parameter y is variationally optimized. As an example, for 
the N3LO-Idaho potential, it can be chosen in the interval 
6-8 fin" 1 . 

The functions J/jc)(i3 (p) ) are written as 



'[Si®; 



1 



[T 2 



1 



J7T- 



(ID 

where the sum is performed over the three even permuta- 
tions. The spin (isospin) of particles i and j are coupled 
to S 2 (Ta), which is itself coupled to the spin (isospin) of 
the third particle to give the state with total spin S (isospin 
T, T z ). The total orbital angular momentum L and the to- 
tal spin S are coupled to the total angular momentum /, J z . 
The functions Y^Z{Qy), having a definite value of L, L z , 
are the HH functions: 



Y (2 {x 2p )®Y h {x lp ) 



N [G] Vp^Wp) 



(12) 

Here Ye^xip) and Y^(x 2p ) are spherical harmonics, N[g\ 
is a normalization factor and <2) ff 1 ,(2 {<p p ) is an hyperspher- 
ical polynomial. The grand angular quantum number G is 
defined as G — 2n + l\ + The notations [G] and {G} of 
Eqs. (O and (Tj) stand for n] and {€ u h, L,S 2 , T 2 , 

S, T; n\, respectively, and fi of Eq. (0 is /j = \G), I. Note 
that each set of quantum numbers {i\,€z, L,S 2 , T 2 , S, T} is 
called "channel", and the antisymmetrization of }f\c\(0 p ' 1 ) 
requires l 2 + 5 2 + T 2 to be odd. In addition, l\ + ( 2 must 
be even (odd) for positive (negative) parity. 

The HH functions having grand angular quantum num- 
ber G constructed in terms of a given set of Jacobi vec- 
tors x\ p , x 2p , defined starting from the particle order i, j, k, 
can always be expressed in terms of the HH functions con- 
structed, for instance, in terms of xi(j,=\),X2(p=i) with the 
same value of G. In fact, the following relation holds 



z 

C'.,C'n' 



11 y^ z (d^ ~> 

u d ,4,n; e[,r 2 ,n' 1 [t [ ,^;n'] ^(p= 1) ' ' 



(13) 



The result is that the antisymmetric functions J/{gj can be 
expressed as a superposition of functions constructed in 
terms of a given order of particles i, j, k, each one having 
the pair i,j in a definite spin and angular momentum state. 
When the two-body potential acts on the pair of particles 

the effect of the projection is easily taken into account. 

The expansion states of Eq. (Q} in momentum space 
can be obtained as follows. Let fik\ p , frk 2p be the conjugate 
Jacobi momenta of the Jacobi vectors, given by 

1 



tik 2p = —(pj-Pi) , 
V2 

[2 1 
p = V 3 (Pt ~ 2 (Pi + ' 



(14) 



Pi being the momentum of the i-th particle. We then define 
a hypermomentum Q and a set of angular-hyperangular 
variables as 



where 



Q= ^k\ p + kl p , 
Qf = [k 2p ,k lp ;ip p ] , 



ki P 

tan <p p = — 

k 2p 



(15) 



(16) 



Then, the momentum-space version of the wave function 
given in Eq. (|9) is 



^ e) ) = 9c,/(0J/(c)(S (G) ), 



(17) 



where J/ |G) (£ (e) ) is the same as J/ {G} (i3 (p) ) of Eq. O with 



x ip -> k ip , and 



g GJ {Q) = (-if \ dp^-J c+2 (Qp)f l {p). (18) 

Jo y 

With the adopted form of flip) given in Eq. (fTUl l. the cor- 
responding functions gc,i(Q) can be easily calculated, and 
they are explicitly given in Ref . ||2TI . 



2.2 The HH Method for Scattering States Below 
Deuteron Breakup Threshold 

We consider here the extension of the HH technique to 
describe N - d scattering states below deuteron breakup 
threshold, when both local and non-local interaction mod- 
els are considered. 

LS J J 

The wave function f ; describing the N — d scatter- 
ing state with incoming orbital angular momentum L and 
channel spin S , parity n = (-) L , and total angular momen- 
tum /, can be written as 



where the sum is restricted to the values £' v {' 2 , and n' such 



that 



\ + t'r, + 2n' = G. The coefficients a ( f\ L ,, ,, , relat- 

1 l t\ ,t2,n; tj ,c 2 ,;r 

ing the two sets of HH functions are known as the Raynal- 
Revai coefficients 11251 . Also the spin-isospin states can be 
recoupled to obtain states where the spin and isospin quan- 
tum numbers are coupled in a given order of the particles. 



ir/LS J J- tr/LS JJ- if/LS JJ- 



(19) 



where T c ~ describes the system in the region where the 

particles are close to each other and their mutual interac- 
ts// 

tions are strong, while T. z describes the relative motion 
between the nucleon and the deuteron in the asymptotic 
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region, where the N - d nuclear interaction is negligible. 

LS J J 

The function W Q % which has to vanish in the limit of 
large intercluster separations, can be expanded on the HH 
basis as it has been done in the case of bound states. There- 
fore, applying Eq. (Q]i, the function W c JJ * can be casted in 
the form 



VP, 



LSJJ 



(20) 



where \¥ p ) is defined in Eqs. (0 and dTTb in coordinate- 
and momentum-space, respectively. 

LS J J 

The function W A z is the appropriate asymptotic solu- 
tion of the relative N - d Schrodinger equation. It is written 
as a linear combination of the following functions, 



^lsjj : ~ ®lsjj ; (p) > 



(21) 



where the sum over p has to be done over the three even 
permutations and 



1=0,2 



[Y,(x 2 p)®S 2 h ® ^ 



®Y L {y„)\ [r 2 ®-] rr . . 



(22) 



Here the spin and isospin quantum numbers of particles 
i and j have been coupled to S 2 and T 2 , with S 2 = 1, 
r 2 = for the deuteron, Wi(X2 P ) is the deuteron wave func- 
tion component in the waves I = 0, 2, y p is the distance 
between N and the center of mass of the deuteron, i.e. 

y p = J^x\ p , Y/(x 2p ) and Yi(y p ) are the standard spher- 
ical harmonic functions, and the functions R^(y p ) are the 
regular (A = R) and irregular (A = I) radial solutions of the 
relative two-body N - d Schrodinger equation without the 
nuclear interaction. These regular and irregular functions, 
denoted as Tii])p) and Ql(\)p) respectively, have the form 



T L (y P ) 



1 



(2L+l)\\q L C L (rj) %p 



g L (y p ) = {2L+\)\\(F l C L m*tep) GL ^ p) , (23) 

where q is the modulus of the N - d relative momentum 
(related to the total kinetic energy in the center of mass 

system by T cm = |j, ju being the N - d reduced mass), 

rj = 2pe 2 /q and £ p = qy p are the usual Coulomb parame- 
ters, and the regular (irregular) Coulomb function Fl(tj, % p ) 
(Gl(j], £ p )) and the factor Cl(tj) are defined in the standard 
way [26 1 . The factor (2L+ 1)! \q L CiSjj) has been introduced 
so that T and Q have a finite limit for q — > 0. The function 
//?({/;?) — [1 _ exp(— 2m/ p )] 2l+1 nas been introduced to reg- 
ularize Gl at small values of y p . The trial parameter b is 
determined by requiring that fuiyp) — > 1 outside the range 
of the nuclear interaction, thus not modifying the asymp- 
totic behaviour of the scattering wave function. A value 
of b — 0.25 far 1 has been found appropriate. The non- 
Coulomb case of Eq. d23l l is obtained in the limit e 1 — > 0. 



In this case, Fi(t],^ p )/^ p and Giij],^ P )l£,p reduce to the 
regular and irregular Riccati-Bessel functions and the fac- 
tor (2L + l)!!Cz,(?7) -> 1 for 77 -> 0. 

LS J J 

With the above definitions, W A z can be written in the 
form 



r 



LSJJ- 



z 

L'S' 



8 W 5ss>&t<S'JJ, + ^LSX'S'^^L'S'J 



(24) 

where the parameters 'R J LSL , s ,{q) give the relative weight 
between the regular and irregular components of the wave 
function. They are closely related to the reactance matrix 
(TC-matrix) elements, which can be written as 

(2L+ 1)!!(2L' + 1)!! q L+L ' +1 C L (? 1 )C L ,(? 1 )'R J LSL , s ,(q)(25) 

By definition of the TC-matrix, its eigenvalues are tan 5ls j, 
6lsj being the phase shifts. The sum over L' and S' in 
Eq. (l24l is over all values compatible with a given J and 
parity n. In particular, the sum over L' is limited to include 
either even or odd values since (-1) L = n. 

The matrix elements ^ sz/s ,(g) and the linear coeffi- 
cients c p occurring in the expansion of ^ SJJz of Eq. (l20b 
are determined applying the Kohn variational principle, 
which states that the functional 



£ = 



2 V3fc 2 



(H-E), 



(26) 



has to be stationary with respect to variations of the trial 

LS J J 

parameters in r N _, z . Here E is the total energy of the sys- 
tem, m is the nucleon mass, and X. is chosen so that 



(27) 



As described in Ref. 0JD, using Eqs. d2Qj and d24]), the 
variation of the diagonal functionals of Eq. d26l l with re- 
spect to the linear parameters c p leads to the following sys- 
tem of linear inhomogeneous equations: 



(28) 



Two different terms D A corresponding to A = R,I are in- 
troduced and are defined as 



D A LS jj^) = W\^lsjj)- 



(29) 



The matrix elements Kj^ T/ S ,(g ) are obtained varying the 
diagonal functionals of Eq. i26\ with respect to them. This 
leads to the following set of algebraic equations 



Zj ^ J LS,L"S"(l) X L'S',L"S" - Y L s,l>s> , 



(30) 



L"S" 

with the coefficients X and Y defined as 

• r LSJJ„I 



xlsz>s> ' £u' rs . u) . 

Ylsms- = -{U% n + Vc JJ " R \-CK,s'jj) ■ (3D 
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Here <p^? JJz,A is the solution of the set of Eq. ( l28b with 
the corresponding D A term. A second order estimate of 
^■ls L'S'W * s gi ven by the quantities [9?^ L , s ,(q)], obtained 
by substituting in Eq. d26l l the first order results. Such second- 
order calculation provides a symmetric reactance matrix. 
This condition is not a priori imposed, and therefore it is a 
useful test of the numerical accuracy. 

In the particular case of q — (zero-energy scattering), 
the scattering can occur only in the channel L = and the 
observables of interest are the scattering lengths. Within 
the present approach, they can be easily obtained from the 
relation 

a Nl j 



V J+ ^. = -KmXjUq). 02) 



An alternative way to solve the scattering problem, used 
when q + 0, is to apply the complex Kohn variational prin- 
ciple to the >S-matrix, as in Ref . |[T8l . 

The approach presented so far for bound and scattering 
states does not have too many differences compared to the 
method presented for instance in Ref. [17|, and known as 
pair-correlated hyperspherical harmonics (PHH) method. 
In fact, in the PHH method a correlation factor is included 
in the HH expansion of Eq. (f20b to take into account the 
strong short-range correlations induced by the realistic two- 
body potentials, like the AVI 8. The presence of correlation 
functions makes the convergence of the expansion much 
faster than in the uncorrected case. However, the PHH 
method cannot be simply implemented when non-local two- 
body interactions are considered, unless the Fourier trans- 
form of the potential is performed. The calculation involv- 

LS J J 

ing W c can be performed with the HH or PHH expan- 
sions in coordinate- or in momentum-space, depending on 
what is more convenient. 



3 Three Nucleon Force Models 

In Ref. [ 13 1 the description of bound states and zero-energy 
states for A — 3, 4 has been reviewed in the context of the 
HH method. In Table [TJ we report results for the triton and 
4 He binding energies as well as for the doublet n-d scatter- 
ing length 2 a nd using the AVI 8 and the N3LO-Idaho NN 
potentials and using the following combinations of two- 
and three-nucleon interactions: AV 1 8+URIX, AV 1 8+TM' , 
N3LO-Idaho+N2LOL and N3LO-Idaho+URIXp. In this 
last model the parameter in front of the spin-isospin in- 
dependent part of the URIX potential has been rescaled by 
a factor of 0.384 to fit the triton binding energy [14] (we 
call this model URIXp). We have considered also the Vuwk 
model, obtained from the AVI 8 interaction with a cutoff 
parameter A = 2.2 fm _1 . The results are compared to the 
experimental values reported in the table. Worthy of notice 
is the recent very accurate datum for 2 a nd l27l . 

From the table we may observe that only the results ob- 
tained using an interaction model that includes a TNF are 
close to the corresponding experimental values. In the case 
of the AVI 8+TM', the strength of the TM' potential has 
been fixed to reproduce the 4 He binding energy and, as can 
be seen from the table, the triton binding energy is under- 
predicted. Conversely, the strength of the URIX potential 



Table 1. The triton and 4 He binding energies B (in MeV), and 
doublet scattering length 2 a„ d (in fm) calculated using the indi- 
cated two- and three-nucleon interactions. The experimental re- 
sults are also reported. 



Potential 


BCH) 


S( 4 He) 




AV18 


7.624 


24.22 


1.258 


N3LO-Idaho 


7.854 


25.38 


1.100 


AV18+TM' 


8.440 


28.31 


0.623 


AV18+URIX 


8.479 


28.48 


0.578 


N3LO-Idaho+N2LOL 


8.474 


28.37 


0.675 


N3LO-Idaho+URIXp 


8.481 


28.53 


0.623 


Vlow-k 


8.477 


29.15 


0.572 


Exp. 


8.482 


28.30 


0.645±0.003±0.007 



has been fixed to reproduce the triton binding energy giv- 
ing too much binding for 4 He. The strength of the N2LOL 
potential has been fixed to reproduce simultaneously the 
triton and the 4 He binding energies whereas the N3LO- 
Idaho+URIXp model overbinds 4 He. These two models 
give a better description of 2 a„d- The Vi ow -k interaction re- 
produces the triton binding energy but overbinds 4 He ap- 
preciably and 2 a nd is not well described. In conclusion a 
simultaneous correct description of the three quantities is 
not achieved by any of the models considered. 

To analyze further this fact, we give a brief description 
of the TM' (or Brazil), URIX and N2LOL models. They 
can be put in the following way: 

W(l,2,3) = aWa(h2,3) + bW b (l,2,3) + dW d (l,2,3) 
W(l,2,3) = aW a (l,2,3) + bW b (l,2,3) + dW d (J,2,3) 
+c D W D (l,2,3) + c E W E (l,2,3). (33) 

Each term corresponds to a different source and has a dif- 
ferent operator structure. The first three terms arise from 
the exchange of two pions between three nucleons. The 
a-term is coming from tiN S -wave scattering whereas the 
fo-term and li-term, which are the most important, come 
from ttN P-wave scattering. The specific form of these 
three terms in configuration space is the following: 

W a (l,2,3) = -^(ti ■ Tz)(cri ■ r 31 )(er 2 • r 23 )y(r 3l )y(r 23 ) 
W/,(l,2,3) = Wq(ti ■ T 2 )[(cri ■ (T 2 )y(r 3l )y(r 23 ) 

+ (ci ■ r 3 i)((T 2 ■ r 23 )(r 31 ■ r 23 )t(r 3 i)t(r 23 ) 
+ (cti • r 3l )(o- 2 ■ r 31 )t(r 3l )y(r 23 ) 
+ (ci ■ r 32 )(a 2 ■ r 32 )y(r 31 )t(r 23 )] (34) 
Wd(l,2,3) = W (t 3 ■ n x t 2 )[(<7 3 ■ <t 2 x cT X )y(r 3x )y{r 23 ) 
+ (o-i ■ r 3 i)((T 2 ■ r 23 )(er 3 ■ r 31 X r 23 )t(r 31 )t(r 23 ) 
+ (ctj • r 3 i)(<T 2 • r 3] x cr 3 )t(r 31 )y(r 23 ) 
+ (C2 ■ r 32 )((T 3 ■ r 32 x (Ti)y(r 3 i)t{r 23 )] , 

with Wo an overall strength. The b- and of-terms are present 
in the three models whereas the a-term is present in the 
TM' and N2LOL and not in URIX. In the first two mod- 
els, the radial functions y(r) and t(r) are obtained from the 
following function 



Mr) 



\2n 1 r 

ml 2n 2 Jo 



dqq' 



Mqr) 

q 2 + ml 



F A (q) (35) 
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where m K is the pion mass and 

y(r) = i/o'W 

(36) 

t(r) = \y'(r) . 

The cutoff function F \ in the TM' or Brazil models is taken 
as [(A 2 -ml)/(A 2 +q 2 )] 2 . In the N2LOL model it is taken as 
exp(-q 4 /A 4 ). The momentum cutoff A is a parameter of the 
model fixing the scale of the problem in momentum space. 
In the N2LOL, it has been taken A = 500 MeV, whereas 
in the TM' model the quantity Ajm n has been varied to 
describe the triton or 4 He binding energy at fixed values of 
the constants a,b and d. In the literature several cases have 
been explored with typical values around A = 5m ff . 

In the URIX model the radial dependence of the b- and 
af-terms is given in terms of the functions 

Y(r) = e-'/xfr 

(37) 

T( r ) = (l + 3/x + 3/x 2 )Y(r){ T 

with x = m n r and the cutoff functions are defined as £y = 
%t = (1 - e~ )> with c — 2.1 firT 2 . This regularization has 
been used in the AVI 8 potential as well. Since the param- 
eters in the URIX model has been determined in conjunc- 
tion with the AVI 8 potential, the use of the same regular- 
ization was a choice of consistency. The relation between 
the functions Y(r), T(r) and those of the previous models 
is 

Y(r) = y(r) + T(r) 

(38) 

T{r) = $t(r) . 

With the definition given in Eq.(l35l). the asymptotic be- 
haviour of the functions /o(r), y(r) and f(r) is: 

3 e~ A 

Mr -> oo) -> — 

x 

3eT x i 1 \ 
i,(r->oo) _-_h + _l (39) 

3 e -*/ 3 3 \ 
t(r -> oo) -> — — i + . 

r z x \ x x l J 

In fact, with the normalization chosen for /o, the functions 
Y and T defined from y and t in Eq. d38l ) and those ones de- 
fined in the URIX model in Eq. d37l > coincide at large sep- 
aration distances. Conversely, they have a different short 
range behavior. 

The last two terms in Eq. d33l correspond to a 2N con- 
tact term with a pion emitted or absorbed (D-term) and to a 
3N contact interaction (/s-term). Their local form, in con- 
figuration space, derived from Ref. (9), are 

W D {\,2,3) = <(T!-T 2 )X 

{(<Tl ■ <x 2 )[i/(r 31 )Z (r 23 ) + Zo(r 31 )y(r 2 3)] 
+ (cri ■ r 31 )(er 2 ■ r 3 i)f(r 3 i)Zo(r23) 
+ (<ti • r 32 )(<T 2 • r 32 )Z (r 31 )f(r 23 )} (40) 
W £ (l,2,3) = W§(n ■ T 2 )Z (r 31 )Z (r 23 ) . 



The constant Wfi, fix the strength of these terms. In 
the case of the URIX model the E'-term is present without 
the isospin operator structure and it has been included as 
purely phenomenological, without justifying its form from 
a particular exchange mechanism. Its radial dependence 
has been taken as Zo(r) = T 2 (r). In the N2LOL model, 
the function Zo(r) is defined as 

Z «W = ^¥A f dqq 2 Mqr)F A {q) (41) 
mi 2n l Jo 

with the same cutoff function used in the definition of /o in 
Eq.([35), F A {q) = exp(-q 4 /A 4 ). In the TM' model the D- 
and fi-terms are absent. 

Each model is now identified from the values assigned 
to the different constants a, b, d, cp, ce- Following Refs. II6I28I . 
in the case of the TM' model, the values of the constants 
have been chosen as a — -0.87 w" 1 , b - -2.58 m n 3 , and 
d = -0.753 TO" 3 ; the strength Wo = (gntx/Smnti) 2 m\ and 
the cutoff has been fixed to A = 4.756 m n in order to de- 
scribe correctly B( 4 He). In Table [T] the calculations have 
been done using these values with g 2 = 197.7, m„ = 139.6 
MeV, m^/m^ = 6.726 (m^ is the nucleon mass) as given in 
the original derivation of the TM potential. As mentioned 
before, this model does not include the D- and E'-terms. 

In the URIX model the b- and cZ-terms are present, 
however with a fix relative value. The strength of these 
terms is: bW = 4 A™ and d = b/4, with A™ = -0.0293 
MeV. The model includes a purely central repulsive term 
introduced to compensate the attraction of the previous 
term, which by itself would produce a large overbinding 
in infinite nuclear matter. It is defined as 

W™ /X (l,2,3) = A R r 2 (r 31 )r 2 (r 23 ) (42) 

with A R = 0.0048 MeV. 

In the N2LOL potential the constants of the a-, b-, d-, 
D- and E'-terms are defined in the following way: 

with a = cim 2 n , b = c 3 /2, d = c 4 /4, and c x = -0.00081 
MeV" 1 , c 3 = -0.0032 MeV" 1 , c 4 = -0.0054 MeV" 1 taken 
from Ref. [16|. The other two constants, = 1.0 and 
ce = -0.029, have been determined in Ref. [9 1 from a fit to 
B( 3 H) and B( 4 He) using the N3LO-Idaho+N2LOL poten- 
tial model. The numerical values of the constant entering 
in W , W° D and W° E are taken as m„ = 138 MeV, F„ = 92.4 
MeV, gA — 1 .29, and the chiral symmetry breaking scale 
A x = 700 MeV. 

In order to analyze the different short range structure of 
the TNF models, in Fig.[T]we compare the non-dimensional 
functions Zo(r), y(r) and T(r) for the three models under 
consideration. In the TM' model using the definition of 
Eq. (l4TT > and using the corresponding cutoff function we can 
define: 



19 International IUPAP Conference on Few-Body Problems in Physics 




This function is showed in the first panel of Fig. Q] as a 
dashed line. From the figure we can see that, in the case of 
the URIX model, the functions Zo(r) and y{r) go to zero as 
r — > 0. This is not the case for the other two models and is 
a consequence of the regularization choice of the Y and T 
functions adopted in the URIX. 

4 Parametrization Study of the Three 
Nucleon Forces 

In this section we study possible variations to the parametriza- 
tion of the TNF models in order to describe the A = 3,4 
binding energies and 2 a„d. 



4.1 Tucson-Melbourne Force 

We first study the TM' potential and we would like to see 
if, using the AV18+TM' interaction, it is possible to repro- 
duce simultaneously the triton binding energy and the dou- 
blet n - d scattering length for some values of the param- 
eters. The a-term gives a very small contribution to these 
quantities, therefore, in the following analysis we maintain 
it fixed at the value a = -0.87 m~ l . In Fig. [2 left panel, the 
doublet n - d scattering length is given as a function of the 
parameter b (in units of its original value b = -2.58 m~ 3 ) 
for different values of the cutoff A (in units of m n ). The 
box in the figure includes values compatible with the ex- 
perimental results. The value of the constant d has been 
fixed to reproduce the triton binding energy. The corre- 
sponding values of the parameter d (in units of its original 
value d = -0.753 m n 3 ) are given in the right panel as a 
function of b. Each point of the curves in both panels cor- 
responds to a set of parameters that, in connection with the 
AV 1 8 potential, reproduces the triton binding energy. The 
variations of the parameters given in Fig.|2]do not exhaust 
all the possibilities. However we can observe that, with the 
AV18+TM' potential, there is a very small region in the 
parameter's phase space available for a simultaneous de- 
scription of the triton binding energy and the doublet scat- 
tering length. This small region corresponds to a big value 
of b and d results to be almost zero. Moreover, the value of 
the cutoff A around 3.8m„- is smaller than the values usually 
used with the TM' potential (A ~ 5m„). 

To be noticed that, for negative values of the parame- 
ters a, b and d, the TM' potential is attractive. It does not 
include explicitly a repulsive term. Added to a specific NN 
potential that underpredicts the three-nucleon binding en- 
ergy, it supplies the extra binding by fixing appropriately 
its strength. However, as mentioned in the Introduction, the 



scattering length is sensitive to the balance between the at- 
tractive part and the repulsive part of the complete interac- 
tion. Therefore, it seems that supplying only an attraction, 
fixed to reproduce the triton binding energy, in the case of 
the TM' interaction it is difficult to reproduce correctly this 
balance. 

As discussed before, the TM' potential is a modifica- 
tion of the original TM potential compatible with chiral 
symmetry. At the same order (next-to-next-to-leading or- 
der) in the chiral effective field theory the D- and E'-terms 
appear (see Ref. J8[ and references therein) as given in 
Eq.(l33l). Here we introduce the following additional term 
to the TM' potential based on a contact term of three nu- 
cleons 

W™(1, 2, 3) = W E J] Z™(r 3l )Z™(r 23 ) . (45) 

eye 

This term is similar to the repulsive term of the URIX 
model and, for the sake of simplicity, we do not include the 
(n • Tz) operator. The function Z™ is a positive function, 
therefore, for positive values of eg, the new term is repul- 
sive. We include it in the following analysis of the TM' 
potential. The analysis of the new term is given in Fig. [3] 
In the left panel the doublet n-d scattering length is given 
as a function of the parameter b (in units of its original 
value b = -2.58 m~ 3 ) for different values of the strength 
of the Wj M -term. The value of the cutoff A has been fixed 
to 4.8 m„. The box in the figure includes values compat- 
ible with the experimental results. Moreover, the value of 
the constant d has been fixed to reproduce the triton bind- 
ing energy. The corresponding values of the 4 He binding 
energy, B( 4 He), is given in the right panel. 

Comparing the left panels in Figs. [2] and [3] the effect 
of the new term is clear. In Fig. [2] we see that using A = 
4.8 m n , 2 a n d is not well reproduced. Conversely, in Fig. [3] 
the inclusion of the new term moves this curve in the cor- 
rect direction and with values of its strength around ce = 
1.6 it is possible to reproduce the experimental value of 
2 a n d. There is also an improvement in the description of 
B( 4 He). In fact, the AV18+TM' model with A = 4.8 m„ 
reproduces the triton binding energy as can be seen from 
Fig. 12 However it predicts B( 4 He) = 28.55 MeV, which is 
slightly too high. With the W™-term, at c E = 1.6, the 
description of B( 4 He) improves. For example with b = 
-3.87 m~ 3 , d = -3.375 m~ 3 and A = 4.8 m n , we ob- 
tain B( 4 He) = 28.36 MeV, very close to the experimental 
value. 



4.2 UrbanalX Force 

In the following we analyze the URIX potential which has 
two parameters, A™ and Ar. In this model the strength 
of the (i-term was related to the strength of the b-Xerm as 
b = Ad. The original values of the parameters were fixed 
in Ref. [4| in conjunction with the AVI 8 NN potential 
and, from Table Q] we observe that the model correctly 
describes the triton binding energy. However, it overesti- 
mates B( 4 He) and underestimates 2 a„d- In order to improve 
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Fig. 1. The Z (r), y(r) and T(r) functions as functions of the interparticle distance r for the URIX (solid line), TM' (dashed line) and 
N2LOL (dotted line) models. 




the description of these quantities, we have varied the con- 
stants A^T, Ar and the relative strength = d/b be- 
tween the b- and li-terms. For a given value of A^Y, the 
values of Ar and has been chosen to reproduce B( 3 H) 
and 2 a n d. The results are given in Fig. [4] In panel (a), A™ 
is given as a function of with Ar varying from 0.0176 
MeV at A™ = -0.02 to 0.0210 MeV at A™ = -0.050 
MeV. These values of Ar are more than three times greater 
than the original value. In panel (b) and (c) the results for 
2 a„d and B( 4 He) are given respectively. The latter has not 
been included in the determination of the parameters, how- 
ever we observe a rather good description in particular for 
values of D™ > 0.7. 



With a modification of the parameters in the URIX 
force, we were able to describe B( 3 H), 2 a„d and B( 4 He). 
This has been achieved with a substantial increase of the 
repulsive term. Also D^Y is quite far from its original value. 
For example, at the original value of A™ = -0.0293 MeV, 
the relative strength is D™ = 1 and A R = 0.0181 MeV. 
This is four times and more than three of the original val- 
ues, respectively. As diminishes, Ar tends to increase 
further with the consequence that the mean value of the 
repulsive part of W results to be more than three times 
the original AV18+URIX value. This is compensated by 
a lower mean value of the kinetic energy. A further analy- 
sis of the effects of the new parametrizations is done in the 
next section studying selected p - d polarization observ- 
ables. 
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Fig. 3. The doublet scattering length £(„_,/ as a function of the parameter b of the TM' potential including the W™-term, for different 
values of the strength c E (right panel). The corresponding values of S( 4 He) (left panel). 
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Fig. 4. (a) The relative strength as a function of A^. In each point of the curve the triton binding energy and 2 a„ d are well described, 
(b) Values of 2 a nd for the seven combinations of the parameters indicated as solid points in panel (a), (c) The corresponding predictions 
for B( 4 He). The crosses indicate the results using the parameters defined in the URIX model 



4.3 N2LOL Force 

The parameters c\, and c<\ of the N2LOL have been 
taken from the chiral N3LO NN force of Ref. [ 16 1, whereas 
the cq and ce parameters have been determined in Ref. (9), 
in conjunction with that NN force, by fitting B( 3 H) and 
B( 4 He). Here we are going to use the N2LOL force in con- 
junction with the AVI 8 NN interaction, so we have to mod- 
ify its parametrization since the amount of attraction to be 
gained is now different (see Table[TJ. Moreover, the modifi- 
cation has to be done in such a way that B( 3 H) and 2 a„d are 



well reproduced. As an example, in Fig. [5] 1 a„d is shown 
as a function of the parameter C3 (in units of its original 
value c 3 = 0.0032 MeV" 1 ) fixing c D = 0.4, c E = 0.1 and 
varying c\ in order to reproduce B( 3 H). With the values 
c 3 = -0.0048 MeV" 1 , c 4 = 0.0043 MeV" 1 , 2 a nd fall inside 
the box and matches the experimental value. In this case, 
the 4 H binding energy results B( 4 H) = 28.36 MeV. 
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5 Polarization observables with the new 
parametrizations 



In the previous section we have analyzed different parametriza- 
tions of the TM\ URIX and N2LOL TNFs determined in 
conjunction with the AVI 8 NN potential. With the new 
parametrizations the three quantities under observation, B( 3 H), 
1 a n d and fi( 4 He), are well reproduced. However, some sub- 
stantial modifications to the first two models were neces- 
sary. In the case of the TM' interaction, we found nec- 
essary to include a repulsive term. In the analysis of the 
URIX interaction, the strength of the repulsive term re- 
sulted to be more than three times larger. In the case of the 
N2LOL interaction, a minor adjustment of the parameters 
was necessary. Now we would like to analyze the effects of 
the new parametrizations in observables that are not cor- 
related to the binding energies or to 2 a ni j. Some polariza- 
tion observables in p - d scattering have this characteris- 
tic, in particular the vector and tensor analyzing powers. 
In Fig. [5] the differential cross section dcr/dQ, the vector 
polarization observables A y and iT\ \ and the tensor polar- 
ization observables 7*20, Ti\ and 7^2 are shown at the lab- 
oratory energy E\ a \, = 3 MeV, for the different potential 
models. As a reference we use the AV18+URIX interac- 
tion given in the figure as a blue line. In the figure, the other 
three curves corresponds to particular parametrizations of 
the models that reproduce 2 a„d and B( 3 H) and approxi- 
mate, as much as possible, B( 4 He). The parametrizations 
of the models selected for the figure are the following: 
the AV18+URIX* model is defined with A™ = -0.0293 
MeV, D™ = 1 and A R = 0.018 MeV In the AV18+TM* 
model we have used a = -0.87 m K l , b - -9.804 m~ 3 , d 

11657 m; 3 , c e = 1, and A = Am„. In the AV18+N2LO* 
model the parametrization corresponds to c\ = -0.00081 
MeV" 1 (its original value), c 3 = -0.0048 MeV" 1 , c 4 = 
-0.0043 MeV" 1 , c D = 0.4 and c E = 0.1. From the figure 
we can observe that the models describe equally well the 
differential cross section and the tensor analyzing powers 
T20, 722- Differences are observed in the vector analyzing 
powers A y and iT\\. Taking as a reference the results of 
the AV18+URIX model, in both cases the AV18+URIX* 
model produces a noticeable worse description whereas 
the AV18+N2LOL* slightly improves the description. The 
new parametrizations of the TNF models overpredict 7^1 
in all cases, in particular the AV18+TM* model. 



6 The Kohn Variational Principle in terms 
of Integral Relations 

Recently two integral relations have been derived from the 
KVP [22|. It has been shown that starting from the KVP, 
the tangent of the phase-shift can be expressed in a form 
of a quotient where both, the numerator and the denomi- 
nator, are given as two integral relations. Let us first con- 
sider a two-body system interacting through a short-range 
potential V(r) at the center of mass energy £ in a rela- 
tive angular momentum state I = 0. The solution of the 
Schrodinger equation in configuration space (to is twice 
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Fig. 5. 1 a„d as a function of the C3 parameter in the N2LOL 
model. 



the reduced mass), 

h 2 



( V z + V -E)¥(r) = 

TO 



(46) 



can be obtained after specifying the corresponding bound- 
ary conditions. For E > 0, with k 2 = E/(h 2 /m) and assum- 
ing a short-range potential V, f (r) = (p(r)l yfAn and 



<p(r — * 00) 



sin(Ar) cos(fcr) 
A h B 



kr 



kr 



(47) 



With the above normalization, the solution ¥ verifies the 
following integral relations: 



- § < ¥\H - E\F >= B with F = 

f<¥\H- E\G >= A with G = 
land = 4 . 




Explicitly they are 



H 2 Vl 



drsm(kr)V(r)[r<p(r)] = B 



to f M 0(0) 

drcos(kr)V(r)[r<f>(r)] + ^ = A , ( 4 9) 

n 2 yfkJo \[k 



where in the last integral we have used the property V 2 (l / r) 
-An5(r). 

In practical cases the solution of the Schrodinger equa- 
tion is obtained numerically. Then, tan 6 is extracted from 



19 International IUPAP Conference on Few-Body Problems in Physics 




<p(r) analyzing its behavior outside the range of the poten- 
tial. The equivalence between the extracted value and that 
one obtained from the integral relations defines the accu- 
racy of the numerical computation. A relative difference 
of the order of 10~ 7 of the two values is usually achieved 
using standard numerical techniques to solve the differen- 
tial equation and to compute the two one-dimensional inte- 
grals. To be noticed the short range character of the integral 
relations. This means that the phase-shift is determined by 
the internal structure of the wave function. 

The last relation in Eq. d49b shows a dependence on 
the value of the wave function at the origin. It could be 
convenient to eliminate this explicit dependence since the 
numerical determination of <p(0) might be problematic, as 
we will show. To this end we introduce a regularized func- 



tion G = f reg G with the property |G(r = 0)| < oo and 
G = G outside the interaction region. A possible choice is 



G= A /±££^2(l- e -r 



(50) 



where the regularization function f reg — (1 -e~ yr ) has been 
introduced with y being a non linear parameter which will 
be discussed below. Values verifying y > 1/ro, with ro 
the range of the potential, could be appropriate. The reg- 
ularized function G (as well as the irregular function G), 
verifies the normalization condition 



^ [< F\H -E\G>-< G\H - E\F >] = 1 



(51) 
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Therefore the second integral relation in Eq. ( l48l l remains 
valid using G in place of G, 



— < <F\H-E\G >= A , 

7i 2 



(52) 



with the explicit form: 



n 2 Vfe 



f 

Jo 



drcos(kr)V(r)[rcf>(r)] + L = A (53) 



where in l y all terms depending on y, introduced by f reg , 
are included. Comparing Eq. (l53l to Eq. d49l we identify 
l y = 0(0)/ jk. 

In the following we demonstrate that the relation tan 6 = 
B/A, which is an exact relation when the exact wave func- 
tion W is used in Eq. (l48l i. can be considered accurate up to 
second order when a trial wave function is used, as it has a 
strict connection with the Kohn variational principle. 

The connection of the integral relations with the KVP 
is straightforward. Defining a trial wave function W t as 



W t =W c +AF + BG 



(54) 



with W c — > as r — > oo, the condition W, — > AF + B G as 
r — > oo is fulfilled. The KVP states that the second order 
estimate for tan 6 is 

[tan<5] 2 "" = tand-^r < (1 /A)Y,\H - E\(l /A)Y, > . (55) 

The above functional is stationary with respect to varia- 
tions on W c and tan<S. Without loosing generality W c can 
be expanded in a (square integrable) complete basis 



(56) 



The variation of the functional with respect to the linear 
parameters a„ and tan 6 leads to the following equations 



<<h\H-E\¥, >=0 
<G\H-E\% >=0 . 



(57) 



To obtain the last equation, the normalization relation of 
Eq. ( BTT l has been used. From these two equations, f c and 
the first order estimate of the phase shift (tand) 1 can be 
determined. To be noticed that the first equation implies 
< W C \H — E\W, >- 0. Furthermore, from the general rela- 
tion (m/h 2 ) [< W,\H -E\G>-< G\H - E\ V, >] = A, and 
using the second equation in Eq. ( |57] |. the following inte- 
gral relation results 



-r < %\H - E\G >= A . 



(58) 



Replacing the two relations of Eq.(l57t into the func- 
tional of Eq.(l55l). a second order estimate of the phase shift 
is obtained 



[tan 6] 1 



(tantf) 1 "-^ <F\H-E\(\IA)W t > . (59) 



Multiplying Eq. 



by A one gets 



B =B 1 - — r < F\H - E\W,> 



(60) 



On the other hand, a first order estimate for the coefficient 
B can be obtained from the general relation 

^ [< F\H-E\W t > - < W,\H-E\F >] = B l " . (61) 
n 

Therefore, replacing Eq.doTTl in Eq.d60l>, a second order in- 
tegral relation for B is obtained. The above results can be 
summarized as follow 

v ,d m 
B =-— < ¥ t \H-E\F> 

n £ 

A = 5 < W,\H-E\G > 

2 nd , 



[tan<5] 2 = B /A 



(62) 



These equations extend the validity of the integral rela- 
tions, given in Eq.(l48l> for the exact wave functions, to trial 
wave functions. To be noticed that F, G are solutions of the 
Schrodinger equation in the asymptotic region, therefore 
(H - E)F — > and (H - E)G — > as the distance between 
the particles increases. As a consequence the decomposi- 
tion of W, in the three terms of Eq. d54l > can be considered 
formal since, due to the short-range character of the rela- 
tion integrals, it is sufficient that the trial wave function be 
a solution of (H — E)W t = in the interaction region, with- 
out an explicit indication of its asymptotic behavior. This 
fact, together with the variational character of the relations 
allows for a number of applications to be discussed in the 
next sections. 



7 Integral Relations for A = 2,3 systems 

Applications of the integral relations to systems with A = 
2, 3 are given. We first consider the following central, s- 
wave gaussian potential 



V(r) = -V exp(-r 2 /r 2 ) 



(63) 



with V = -51.5 MeV, r = 1.6 fm and h 2 /m = 41.4696 
MeV fm 2 . This potential has a shallow L = bound state 
with energy E 2B = -0.397743 MeV. 

In the A = 2 system, the orthogonal basis 



4 2) (z)exp-(z/2) 



(64) 



with X m a (normalized) Laguerre polynomial and z = Br, 
being B a nonlinear parameter, is used to expand the wave 
function of the system 



z 



(65) 



The Schrodinger equation is transformed to an eigenvalue 
problem that can be solved for different values of the di- 
mension M of the basis. The variational principle states 
that 

E = (%\H\W ) > E 2B , (66) 
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with the equality obtained for M — > oo. The nonlinear pa- 
rameter B can be fixed to make improve the convergence 
properties of the basis. In fact, for each value of M there is 
a value of B that minimizes the energy. Increasing M, the 
minimum of the energy becomes less dependent on B re- 
sulting in a plateau. Increasing further the dimension of the 
basis, the extension of the plateau increases as well, with- 
out any appreciable improvement in the eigenvalue, indi- 
cating that the convergence has been reached up to certain 
accuracy. At each step Wq represents a first order estimate 
of the bound state exact wave function. 

In the proposed example the system has only one bound 
state. So, with proper values of M and the diagonaliza- 
tion of H results in one negative eigenvalue Eq and M-1 
positive eigenvalues Ej (j = 1, M— 1). The correspond- 
ing wave functions 



M-1 

z< 



Table 2. The two-nucleon bound state Eq and the first three pos- 
itive eigenvalues Ej (j = 1,3), as a function of the number of 
Laguerre polynomials M. The second order estimates, [tan<^] 2 , 
obtained applying the integral relations are given in each case and 
compared to exact results, tan<5j. 



M 


10 


20 


30 


40 


Eq 


-0.395079 


-0.397740 


-0.397743 


-0.397743 


Ei 


0.536349 


0.116356 


0.048091 


0.026008 


[tancS,] 2 "" 


-1.507280 


-0.622242 


-0.392005 


-0.286479 


tan<5i 


-1.522377 


-0.621938 


-0.392021 


-0.286480 


E 2 


1.984580 


0.449655 


0.190019 


0.103503 


[tan(5 2 ] 2 "' 


-5.919685 


-1.353736 


-0.812313 


-0.584389 


tan<S 2 


-5.703495 


-1.354691 


-0.812270 


-0.584388 


E 3 


4.512635 


0.994433 


0.423117 


0.231645 


[tancy 2 "" 


13.998124 


-2.451174 


-1.302799 


-0.908128 


tan S3 


12.684474 


-2.448343 


-1.302887 


-0.908131 



,M- 1 



(67) 



are approximate solutions of (H - Ej) Wj = in the inter- 
action region. As r — > 00 they go to zero exponentially and 
therefore they do not represent a physical scattering state. 
The negative energy £0 and the first three positive energy 
eigenvalues {Ej, j — 1, 3) are shown in Fig.|7]as a function 
of f3 in the case of M — 40. We observe the plateau already 
reached by £0 for the values of B showed in the figure. 
We observe also the monotonic behavior of the positive 
eigenvalues toward zero as B decreases. The corresponding 
eigenvectors can be used to compute the integral relations 
of Eq. ( 1621 and to calculate the second order estimate of 
the phase-shifts 8j at the specific energies Ej. This analy- 
sis is shown in Table [2]in which the non linear parameter B 
of the Laguerre basis has been fixed to 1 .2 frrT 1 . In the first 
row of the table the ground state energy is given for differ- 
ent values of the number M of Laguerre polynomials. The 
stability of Eq at the level of 1 keV is achieved already with 
M = 20. For a given value of M, Ej, with j - 1, 2, 3, are 
the first three positive eigenvalues. The eigenvectors cor- 
responding to positive energies approximate the scattering 
states at the specific energies. Since the lowest scattering 
state appears at zero energy, none of the positive eigenval- 
ues can reach this value for any finite values of M. Defining 
k 2 j = fpEj, the second order estimate for the phase shift at 
each energy and at each value of M is obtained as 



' kj sin(feyr) 




%<¥j\H-E\Fj>=BjWiih Fj 



§ < Tj\H - E\Gj >= Aj with Gj = f t 

[tan 5 jf d = Bj/Aj. 

On the other hand, as we are considering the A = 2 
system, at each specified energy Ej the phase shift tan Sj 
can be obtained by solving the Schrodinger equation nu- 
merically. The two values, [tancSj] 2 and tanSj, are given 
in the Table [2] at the corresponding energies as a function 
of M. We observe that, as M increases, the relative differ- 
ence between the variational estimate and the exact value 
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— E =-0.397743 MeV 
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Fig. 7. The two-nucleon bound state energy Eq and the first three 
positive eigenvalues Ej as a function of /? in the case of M = 40 



reduces, for example at M — 40 is around 10 -6 . In fact, 
as M increases, each eigenvector gives a better representa- 
tion of the exact wave function in the internal region and 
the second order estimates, [tiandj] 2 approach the exact 
result. 

In a different application, the integral relations can be 
used to calculate the phase-shift of a process in which the 
two particles interact through a short range potential plus 
the Coulomb potential, imposing free asymptotic condi- 
tions to the wave function. As an example we use the same 
two body potential used in the previous analysis and add 
the Coulomb potential: 



V(r) = -V exp-(r/r r + 



(69) 



For positive energies and I = 0, the wave function behaves 
asymptotically as 



^)( r _> oo) = AF c (r) + BG c (r) 



(70) 
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with F c (r), G c (r) the regular and irregular Coulomb func- 
tions, respectively. The phase-shift is tanc) c = B/A. The 
KVP remains valid when the long range Coulomb potential 
is considered and its form in terms of the integral relations 
results: 

- § < Wl c) \H - E\F C >= B 
§ < ^ c) \H - E\G C >=A 

[tan<y 2 "" = f . (71) 

with G c = frecjGc and *Pf' a trial wave function behav- 
ing asymptotically as ¥' (c) . Since (H - E)\F C > and (H — 
E)\G C > go to zero outside the range of the short range po- 
tential, the integrals in Eq. dTTb are negligible outside that 
region. Therefore, for the computation of the phase-shift 

(c) (c) 

it is enough to require that iPJ verifies (H - E) T) = 0, 
inside that region. To exploit this fact, we introduce the 
following screened potential: 

2 

V sc (r) = -V exp [-(r/r ) 2 ] + [e^'"'"] - . (72) 

For specific values of n and r sc it has the property of being 
extremely close to the potential V(r) of Eq. (|69| i for r < ro, 
with ro the range of the short range potential. The screen- 
ing factor e~ (f /' v) " cuts the Coulomb potential for r > r sc . 
Using the potential V sc to describe a scattering process, the 
wave function behaves asymptotically as 

YnsJr -» oo) = AF(r) + BG{r) (73) 

with F, G from Eq. (1681 1. since V sc is a short range potential. 
Solving the Schrodinger equation for this potential, it is 
possible to obtain the wave function W n<rx for different val- 
ues of n and r sc . This wave function can be considered as a 
trial wave function for the problem in which the Coulomb 
potential is unscreened. Accordingly it can be used as in- 
put in Eq. (TTTb to obtain a second order estimate of the 
Coulomb phase-shift, 

- f < 'Fn,rJH-E\F c >=B 

§<¥ n>r JH-E\G c >=A 

[tancy 2 "" = f (74) 

where in H the unscreened Coulomb potential is consid- 
ered. This estimate depends on n and r JC as the wave func- 
tion does. In Fig.|8]the second order estimate [tanc) e ] 2 "' is 
shown as a function of r sc for different values of n. The 
straight line is the exact value of tani5 f obtained solving 
the Schrodinger equation. We can observe that for n > 4 
and r sc > 30 fm the second order estimate coincides with 
the exact results. In this example the integral relations de- 
rived from the Kohn Variational Principle have been used 
to extract a phase-shift in presence of the Coulomb poten- 
tial using wave functions with free asymptotic conditions. 

Finally an application of the integral relations to the 
A = 3 system is discussed. To this end we give the gen- 
eralization of the integral relations to the case in which 




Fig. 8. The two-nucleon second order estimate [tan5 c ] 2 "' as a 
function of r sc for different values of n. As a reference the exact 
value for tan d c is given as a straight line. 



more than one channel is open. The coefficients A and B 
of Eq. (l68l i correspond to matrices 

B iJ = ~<'F i \H-E\Fj> 

A tJ = ^< nH - E\Gj > 

R 2 " d =A- l B. (75) 

with R 1 '"' the second order estimate of the scattering ma- 
trix whose eigenvalues are the phase shifts and the indices 
(z, j) indicate the different asymptotic configurations acces- 
sible at the specific energy under consideration. We con- 
sider p - d scattering at E\ a b - 3 MeV using the AVI 8 
potential in the J = l/2 + state. The corresponding scat- 
tering matrix is a 2 x 2 matrix. The corresponding phase- 
shift and mixing parameters have been calculated using the 
PHH expansion and are given in Table [3] From the previ- 
ous discussion we have shown that it is possible to solve an 
equivalent problem with a screened Coulomb potential, so 
with free asymptotic conditions, and then use the integral 
relations to extract the scattering matrix corresponding to 
the unscreened problem. This has been done using Eq. ( |75| ) 
and the results are given in Table [3] using r sc = 50 fm and 
n sc = 5. We observe a complete agreement between the 
two procedures. 



8 Conclusions 

Stimulated by the fact that the commonly used TNF mod- 
els do not reproduce simultaneously the triton and 4 He 
binding energy and the n - d doublet scattering length, we 
have analyzed possible modifications of some of the TNF 
models usually used in the description of light nuclei: the 
TM' and the URIX models. We have also considered the 
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Table 3. Phase-shift and mixing parameters for p — d scattering at 
E, ab = 3 MeV using the AVI 8 potential. Results using the PHH 
expansion (second column) and using the integral relations (last 
column) 





p - d 


Int.Rel. 


4 £>l/2 


-3.563° 


-3.562° 


2 S 1/2 


-32.12° 


-32.12° 


771/2+ 


1.100° 


1.101° 



recent N2LOL model. In each of these models we have var- 
ied the original parameters so as to improve the description 
of the mentioned quantities. Furthermore we have studied 
the description of some p - d polarization observables at 
Eiab = 3 MeV. We have observed that the modification of 
the URIX produces a worse description of the vector po- 
larization observables due to the artificial increase of the 
strength of the repulsive term. The analysis of the TM' 
model has put in evidence the necessity of including a re- 
pulsive term. In the case of the N2LOL model a fine tuning 
of the parameters was possible in order to have an accept- 
able description of the triton and 4 He binding energies and 
the n - d doublet scattering length. Moreover, in the po- 
larization observables we observe an improvement in the 
vector analyzing powers and a slightly worse description 
of 72i. From this analysis we have established a connec- 
tion between the short-range structure of the TNF and the 
polarization observables at low energies. 

In a different application, we have discussed the use of 
the integral relations derived from the KVP in the descrip- 
tion of scattering states. Firstly we have shown the use of 
bound state like wave functions to compute the scattering 
matrix and, in the case of charged particles, the possibility 
of computing phase-shifts using scattering wave functions 
with free asymptotic conditions, obtained after screening 
the Coulomb interaction. Both problems are of interest in 
the study of light nuclei. 
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